Spin currents and magnetoresistance of graphene-based magnetic junctions 
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Using the tight-binding approximation and the nonequilibrium Green's function approach, we 
investigate the coherent spin-dependent transport in planar magnetic junctions consisting of two 
ferromagnetic (FM) electrodes separated by a graphene flake (GF) with zigzag or armchair interfaces. 
It is found that the electron conduction strongly depends on the geometry of contact between the 
GF and the FM electrodes. In the case of zigzag interfaces, the junction demonstrates a spin-valve 
effect with high magnetoresistance (MR) ratios and shows negative differential resistance features 
for a single spin channel at positive gate voltage. In the case of armchair interfaces, the current- 
voltage characteristics behave linearly at low bias voltages and hence, both spin channels are in on 
state with low MR ratios. 
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I. INTRODUCTION 

Spin-dependent transport through magnetic nanoj unc- 
tions is currently receiving increased attention owing to 
the expectation that nanoscale size constraints may max- 
imize the potentially useful magnetoresistance effects. 
Among various types of nanoscale devices, carbon-based 
nanostructures are attractive for spin-polarized electron- 
ics [l|, because the relatively weakness of spin-orbit and 
hyperfme interactions should lead to long spin coherence 
times Q. In particular, the spin-dependent transport 
studies on carbon nanotubes coupled to ferromagnetic 
Co electrodes have shown that nanotubes exhibit a con- 
siderable tunnel MR effect [1, @ . A Cgo molecule is also 
an interesting nanocarbon unit, and recently the spin- 
dependent electron transport for the granular film con- 
sisting of Co nanoparticles embedded in a Co-Ceo com- 
pound matrix has been studied, and the maximum MR 
ratio of about 30% at low bias voltages was reported 

Another related carbon-based structure is graphene 
which demonstrates many useful electronic properties: 
high carrier mobility, long spin relaxation times and 
lengths, extreme flexibility and stability, and metallic 
properties controllable by a gate electrode [U [T3| • More- 
over, due to the flat structure, graphene seems to be eas- 
ier to manipulate than carbon nanotubes. Furthermore, 
interesting properties appear if the graphene is patterned 
into ribbon like geometry called the graphene nanoribbon 
(GNR). For instance, they can be either semiconducting 
with a size dependent gap or metallic [TTj ] . Transport 
properties of GNRs are expected to depend strongly on 
whether they have an armchair or zigzag edge [12j |. In 
GNRs with zigzag edges, transport is dominated by edge 
states which have been observed in scanning tunneling 
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microscopy [TJ]. These states are expected to be spin- 
polarized and make zigzag GNRs attractive for spintron- 

icsiaQ. 

Based on first-principles calculations, Son et al. fl5| 
predicted that the zigzag GNRs become half-metallic 
when an external transverse electric field is applied. Us- 
ing such calculations, Kim and Kim [l6| showed that 
zigzag GNRs connected to two ferromagnetic (FM) elec- 
trodes exhibit very large values of magnetoresistance. 
Using the tight-binding approximation, Bery and Fertig 
[Ti| studied the magnetoresistance of GNR-based spin 
valves in the infinite width limit. They reported a feeble 
magnetoresistance due to the weak dependence of the 
graphene conductivity on the electronic parameters of 
the FM leads. The other spin-dependent properties of 
graphene such as spin field effect transistor [18} , spin Hall 
effects [Til [2(J, and magnetic ordering in graphene [2lJ 
have also been predicted. 

Experimentally, based on the nonlocal magnetoresis- 
tance measurements, spin injection into a graphene thin 
film connected to FM electrodes has been successfully 
demonstrated [12, [H, [24j] and the possibility of spin 
transport and spin precession over micrometer-scale dis- 
tances at room temperature reported [22] . Hill et al. 
[25} have observed large magnetoresistances of several 
hundred ohms in a spin-valve device where a 200 nm 
wide graphene wire is contacted by two soft magnetic 
NiFe electrodes. Also, Wang et al. |26| have investigated 
magnetotransport properties of spin valves consisting of 
graphite flakes connected to FM electrodes, and the mag- 
netoresistance values up to 12% at 7 K were observed 
when an ultrathin MgO tunnel barrier was inserted at 
the FM/graphite interface. These experiments clearly 
demonstrate the possibility of spintronics applications for 
graphene. 

In this paper, we present a theory to investigate the 
spin-dependent properties of planar magnetic junctions 
consisting of a graphene flake (GF) connected by two 
FM electrodes in the presence of bias and gate volt- 
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FIG. 1: Schematic view of two FM/GF/FM magnetic junc- 
tions with (a) armchair and (b) zigzag interfaces, at zero gate 
and bias voltages. In both cases, the magnetic alignments of 
the electrodes are shown in the parallel configuration. For 
the antiparallel configuration, the magnetic alignment of the 
right (drain) electrode should be reversed. 

ages. A GF is a finite rectangular GNR with finite 
number of carbon atoms and, therefore, it can be cou- 
pled to the two square-lattice FM electrodes through its 
armchair or zigzag edges [13, Such configurations 

have been shown in Fig. 1. The present approach indi- 
cates that the spin-polarized transport in such magnetic 
junctions is mainly controlled by the molecular field of 
the FM electrodes, the geometry of FM/GF interface, 
bias and gate voltages, and the strength of FM/GF cou- 
pling. The paper is organized as follows. In Sec. II, 
we present the model and formalism for spin transport 
through FM/GF/FM magnetic junctions where the in- 
terfaces between FM leads and GF can be armchair or 
zigzag. Numerical results and discussions for transmis- 
sion coefficients, spin currents and magnetoresistance are 
presented in Sec. III. We conclude our findings in Sec. 
IV. 



II. MODEL AND FORMALISM 

We consider a system which consists of a rectangu- 
lar GF with armchair and zigzag edges, and two square- 
lattice FM electrodes. The GF can be coupled to the 
magnetic electrodes in such a way that their interface is 
armchair or zigzag. The two models for such structures 



are shown schematically in Fig. 1. Since the electron 
conduction is mainly determined by the graphitic region, 
the electronic structure of this part should be resolved 
in detail. Hence, we decompose the total Hamiltonian of 
the system as 

H = H l + Hc + Hr + Ht , (1) 

where Hl an d Hr are the Hamiltonians of the left (L) 
and right (R) electrodes, He describes the Hamiltonian 
of the graphitic region (GF), and Ht is the coupling of 
the GF to the electrodes. All terms of the total Hamilto- 
nian are described within the single-band tight-binding 
approximation and can be written as 

Hoc = [( £ « ~ ° ' m '")i l i» ~ Kdcc)] Co-Cj'c.a , 

(2) 

He = J]] ( e c$ic,jc ~ ticjc) dl c a dj c ,cr , (3) 

<ic<3c><" 

H T = - S tl ^ c(6 L^ £ W + H - c -) , (4) 

a={L,R} i a ,jc,<? 

where c| a {ci a ,a) and d\ c a (di Cj(T ) create (annihilate) 
an electron with spin a on site i in the electrodes and 
the graphitic region, ti a j a = t^, ti 0t j — tc, and 
ti a ,jc = ^lc if i a nd j are nearest neighbors and zero 
otherwise. e a , the on-site energy of the electrodes, is 
equal to 2tL which acts as a shift in energy, and e c , the 
on-site energy of the GF, is zero except in the presence of 
gate voltage Vg, that shifts the energy levels of the GF. 
Here, —a ■ m a is the internal exchange energy with m Q 
denoting the local magnetization at site i a , and a being 
the conventional Pauli spin operator. 

In this study we set th = tc — t and assume that 
the transport is ballistic [2{|; hence, we set t^c — 0.8 1, 
because the value of t^c should not be smaller than 
the order of t. Such a ballistic approximation is valid 
when the mean free path of the carriers is greater than 
the sample dimensions. High mobility of charge carri- 
ers in graphene, even at room temperature, implies that 
carriers can move long distances without scattering 0]. 
Hence, we expect that the ballistic approximation to be 
appropriate in studying transport in GFs. In addition, 
we assume that the spin direction of the electron is con- 
served in the tunneling process through the graphitic re- 
gion. Therefore, there is no spin-flip scattering and the 
spin-dependent transport can be decoupled into two spin 
currents: one for spin-up and the other for spin-down. 
This assumption is well-justified since the spin diffusion 
length in organics is about 4 nm Q and especially in car- 
bon nanotubes is at least 130 nm Q and in graphene is 
1.5 fim. [3(|, which are greater than the length scale of 
our device. Since the total Hamiltonian does not con- 
tain inelastic scatterings, the spin currents for a constant 
bias voltage, V a , are calculated by the Landauer-Buttiker 
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formula based on the nonequilibrium Green's function 
method [HI]: 

e r°° 

h{V a ) = -J T a (e,V a )[f(e-n L )-f(e-HR)]de, (5) 

where / is the Fermi-Dirac distribution function, (J,l,r — 
Ep ± \eV a are the chemical potentials of the electrodes, 

and T a (e,V a ) = Tr[r L^GaTR^Gl] is the spin-, energy- 
and voltage-dependent transmission function. The spin- 
dependent Green's function of the graphitic region cou- 
pled to the two FM electrodes (source and drain) in the 
presence of the bias voltage is given as 

G a (e,V a ) = [ei~Hc-t L , a (e-eV a /2)-t R , !T (e+eV a l2)]- 1 

(6) 

where Sl i(T and S_r j(T describe the self-energy matrices 
which contain the information of the electronic structure 
of the FM electrodes and their coupling to the GF. These 
can be expressed as T, a . a {e) = Tc. a g a ,a{f)T a ,c where f 
is the hopping matrix that couples the graphene to the 
leads and depends on the geometry of the FM/GF inter- 
face. g a _ a are the surface Green's functions of the uncou- 
pled leads i.e., the left and right semi-infinite magnetic 
electrodes, and their matrix elements are given by 



f-f e + i5- e c 



a ■ m Q + e(k x ,l y ) 



(7) 



where 5 is a positive infinitesimal, s(k x ,l y ) 



2t[cos(k x a) + cos( JV t, ^ 1 )], and 



^/N X (Ny + l 



■. sin(k x Xi) sin( 



N v 



1 



(8) 

Here, l y (= l,...,N y ) is an integer, k x £ [— — , — ], and 
Np with /3 = x, y is the number of lattice sites in the j3 
direction. Using E a ,j, the coupling matrices r Qj(T , also 
known as the broadening functions, can be expressed as 

T a . a = -2Im(S ai(T ). 

When the GF is brought close to an electrode, the 
bonding between them will depend on the arrangement 
and the nature of the atoms at the FM/GF interface. In 
Fig. 1(a), the graphitic region consisting of 5.5 unit cells 
is matched to the FM leads along its armchair edges. In 
this case, the lattice constant of the leads, at, is equal to 
the lattice constant of the GF, ac- On the other hand, 
if two FM electrodes are coupled to the graphene along 
its zigzag edges, we obtain Fig. 1(b) where or, = Viae 
and the GF consists of 3 unit cells. In both junctions, 
there is the same number of carbon atoms at the graphitic 
part, while the number of contact points at the FM/GF 
interfaces has decreased from six (in Fig. 1(a)) to five 
(in Fig. 1(b)). Therefore, one can expect different spin 
transport through the GF, which will be discussed in 
the next section. We should note that, the core of the 
problem lies in the calculation of the spin-dependent self- 
energies Sl.o- and £fl )0 -. In this regard, for a GF with 



N carbon atoms and in the case of contact through n a 
(n z ) carbon atoms at the armchair (zigzag) interfaces, 
only n„ (nj,) elements of the self-energy matrices will be 
non-zero. 



III. RESULTS AND DISCUSSION 

We now use the method described above to study 
the coherent spin-dependent transport and magnetore- 
sistance effect of FM/GF/FM junctions with armchair 
and zigzag interfaces. We have done the numerical cal- 
culations for the case that the direction of magnetiza- 
tion in the left FM electrode is fixed in the +y direction, 
while the magnetization in the right electrode is free to 
be flipped into either the +y or -y direction. We choose 
a GF with iV=66 carbon atoms and set ti, — tc=l eV, 
|m Q |=0.75 eV, E F =0.0 eV, and T=300 K in the calcu- 
lations. 

In Figs. 2 and 3 we depict the spin-dependent trans- 
mission coefficients for the junctions with zigzag (n z =5) 
and armchair (n a =6) interfaces, respectively. Only the 
energy window around the Fermi level has been shown. 
In these figures, the transmission characteristics clearly 
demonstrate the dependence of electron conduction on 
the geometry of contact between the FM electrodes and 
the GF. In the case of parallel alignment of the junction 
with zigzag interfaces, there are non-zero transmissions 
only for spin-up electrons (on state) and the spin-down 
ones are in off state. In the case of antiparallel alignment, 
both spin channels are in off state. This is the so-called 
spin-valve effect or magnetoresistive effect and indicates 
that, the junction with zigzag interfaces can work as a 
switch, or a bit with on and off states. The spin-valve 
devices are promising candidates for systems that may 
transform spin information into electrical signals. 

From the spin-dependent currents of the junction with 
zigzag interfaces it is clear that at low applied voltages 
{Va < 0.5 V) only the spin-up current in the parallel 
alignment is in on state (see Figs. 4(a)-(c)). The current- 
voltage characteristics of such electrons are linear at low 
voltages and display different behaviors at higher volt- 
ages, depending on the value of gate voltage. One in- 
teresting feature of the I-V characteristics is that the 
spin-up current for Vg=0.4 V shows a negative differen- 
tial resistance region which is visible as a decrease in the 
current upon increasing V a [3l( . This behavior is due to a 
gradual disappearance of a resonant level within the volt- 
age window in that region. In addition, the appearance 
of zero current and the necessity of threshold voltage to 
generate finite current flow through the junction in some 
of the I-V curves, arising from the energy mismatch be- 
tween the Fermi energy of the FM electrodes and the 
lowest unoccupied levels of the GF. 

Interestingly, as shown in Fig. 3, the calculated trans- 
mission coefficients of two spin channels in both magnetic 
alignments are significantly larger than the transmission 
spectra of Fig. 2 near the Fermi energy and thus ac- 
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FIG. 2: Transmission coefficients at V" a =0.1 V for 
FM/GF/FM junction with zigzag interfaces in the parallel 
and antiparallel magnetized states. Dashed, solid, and dotted 
curves correspond to Vg = —0.4, 0.0 and 0.4 V, respectively. 



FIG. 3: Transmission coefficients at V" a =0.1 V for 
FM/GF/FM junction with armchair interfaces in the parallel 
and antiparallel magnetized states. Dashed, solid, and dotted 
curves correspond to Vg = —0.4, 0.0 and 0.4 V, respectively. 



cording to Figs. 5(a)-(c) our theory predicts a linear 
behavior in the I-V characteristics at low bias voltages 
for the junction with armchair interfaces. In contrast to 
the spin currents of the junction with zigzag interfaces, 
none of the spin currents exhibits a gap in the I-V char- 
acteristics of Fig. 5, which is due to the existence of 
the resonant states near the Fermi energy in the junction 
with armchair interfaces. 

The main factor of difference in the calculated trans- 
mission spectra arises due to the quantum interference 
effects. In the case of armchair interfaces, the current 
flow mechanism corresponds to the resonant tunneling 
which can be attributed to Fabry-Perot interference of 
electronic states partially reflected from the electrodes, 
due to the finite length of the graphene sample [28[ . In 
the case of zigzag interface, however, the interference be- 
tween electronic waves scattered by the carbon atoms 
along the FM/GF interfaces becomes important, some 
resonances might completely disappear, and the trans- 
mission curve changes (compare Fig. 2 with Fig. 3). 
Therefore, the current flow mechanism is tunneling. The 
physical meaning of the interference effect is that the 
electron waves in the GF which come from the FM elec- 
trode may suffer a phase shift. Thus, a constructive or 
destructive interference in the propagation process of the 
electron through the GF can occur. 

Another interesting feature of the junctions is the mag- 
netoresistance ratio which can be obtained from the I-V 
curves using the definition: MR = (I v — I a )/I p , where 
Ip,a are the total currents in the parallel and antiparallel 
alignments of magnetizations in the FM electrodes, re- 



spectively. In Figs. 4(d) and 5(d) we have shown the 
MR ratios for the two junction at positive and negative 
gate voltages. It is clear that the MR in the junction with 
zigzag interfaces has high values in comparison with the 
other junction. The origin of magnetoresistance effect is 
the difference in the total currents which is related to 
the asymmetry of surface density of states of the FM 
electrodes for spin-up and spin-down electrons and the 
quantum tunneling phenomenon through the GF. In the 
parallel alignment, minority (majority) electrons in the 
left FM electrode go into the minority (majority) states 
of the right one by tunneling through the GF. If, however, 
the two FM electrodes are magnetized in the opposite di- 
rections, the minority (majority) electrons from the left 
electrode seek empty majority (minority) states in the 
right electrode. Consequently, the parallel arrangement 
gives much higher total current through the GF than 
does the antiparallel arrangement in the selected volt- 
age interval except for zero gate voltage at high (low) 
applied voltages for the junction with zigzag (armchair) 
interfaces. 

We also studied the effect of tLC on the spin currents 
and MR ratio. The results showed that in both junc- 
tions, the values of currents increase with increasing t^c 
(from 0.5 to 1 eV), while the MR values do not change 
considerably due to its relation to relative change of total 
currents between both magnetic alignments. 

Now we present a discussion on the spin states at 
the edges of GF. As we pointed out, zigzag edge GNRs 
present half-metallic behavior when homogeneous elec- 
tric fields are applied across the ribbons (in our model, 
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FIG. 4: (a)-(c) Spin currents as a function of applied voltage 
for FM/GF/FM junction with zigzag interfaces. Square and 
circle (up-triangle and down-triangle) curves correspond to 
the spin-up and spin-down currents in the parallel (antipar- 
allel) alignment, (d) MR as a function of applied voltage at 
different gate voltages. 



FIG. 5: (a)-(c) Spin currents as a function of applied voltage 
for FM/GF/FM junction with armchair interfaces. Square 
and circle (up-triangle and down-triangle) curves correspond 
to the spin-up and spin-down currents in the parallel (antipar- 
allel) alignment, (d) MR as a function of applied voltage at 
different gate voltages. 



along the y direction of Fig. 1(a)). This behavior is re- 
lated to the existence of electronic states at the Fermi en- 
ergy spatially localized at the zigzag edges. The ground- 
state spin configuration at zero electric field corresponds 
to an opposite spin orientation across the ribbon between 
the spin-polarized edge states. That is, the spins on the 
outermost sites of each zigzag edge are parallel, resulting 
in one edge polarized with spin up and the other polar- 
ized with spin down, and the total zigzag-GNR spin is 
equal to zero [l5|. In contrast, the armchair edge GNRs 
do not show such spin states due to the delocalized nature 
of the frontier orbitals [32|. Theoretical studies have re- 
vealed that in zigzag GNRs; (i) the critical electric field to 
achieve half-metallicity decreases with the increase of rib- 
bon width, (ii) the ribbon remains half-metallic only at a 
limited range of electric field, and (iii) the half-metallicity 
will be destroyed by a too stron g el ectric field and the 
ribbon becomes nonmagnetic flU l33l |34| . 

From a theoretical point of view, in order to observe 
the spin states due to the edge-localized states in the 
zigzag GNRs, one can use of first-principles calculations, 
which can predict the electrical and magnetic properties 
of a material from the atomic number and mass of its con- 
stituent atoms 0, [SUSHI, HII, or include the electron- 
electron interaction in the total Hamiltonian [3(| \3l\ . 
The structure used in this study is analogous to that 
of Ref. [16]. However, we have used of the single-band 
tight-binding Hamiltonian in the absence of interactions, 
because in this work we are interested in the effects of the 
FM/GF interfaces on the spin transport. Hence, the re- 
sults do not yield the spin states at the edges of zigzag rib- 



bon [Fig. 1(a)]. According to the above discussion, such 
a restriction does not invalidate our model for the zigzag 
edge junction, since the electronic states corresponding 
to the edge states, as it is clear from Fig. 3, are far from 
the Fermi energy. The present theory demonstrates well 
the effect of magnetic states of the electrodes on the elec- 
tron conduction through the GF, and also reveals well the 
role of contact geometry in the spin-polarized transport 
through the two above junctions. We should note that, if 
we include the interactions in our Hamiltonian, the mag- 
netic states corresponding to the edges may affect the 
spin currents and the MR. 

IV. CONCLUSION 

In this work, spin polarized transport of armchair 
and zigzag interfaces between a GF and FM electrodes 
(Fig. 1) is investigated in the tight-binding approxi- 
mation. The results indicate that the atomic arrange- 
ment at the interface between the FM electrode and the 
GF has a dominant effect on the electron conduction 
through the graphene-based magnetic junctions. The 
planar FM/GF/FM junction with zigzag interfaces ex- 
hibits a spin-valve effect with MR values as high as 95% 
and negative differential resistance features for a single 
spin channel. In the case of armchair interface, both spin 
channels arc in on state with low MR ratios. In addition, 
the results of bias and gate voltage dependence of the 
spin currents and MR have been reported. 

In this study, we ignored the effects of spin states at 
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the edges of GF and the magnetic anisotropy of the FM 
electrodes due to the reduced dimensionality and the ge- 
ometry of the system. Although, these factors can affect 
the spin-dependent transport, the present study advances 
the fundamental understanding of intcrfacial effects in 



the graphene-based magnetic junctions and suggests that 
the junction with zigzag interfaces is an interesting can- 
didate for application in the magnetic memory cells and 
spintronic devices. 
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